%%

s = tf('s')

A=1/10;
wd = 10;
w0 = 50;

num = s^2+2*A*wd*s+w0^2
den = s^2+2*wd*s+w0^2

Hs = num/den

bode(Hs)